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Abstract 



Evolution of the z-component of a single spin in the finite cyclic XY spin 1/2 
^ chain is studied. Initially one selected spin is polarized while other spins are com- 

^ pletely unpolarized and uncorrelated. Polarization of the selected spin as a function 

of time is proportional to the autocorrelation function gQ^{t) at infinite temperature. 
The initial perturbation gives rise to two wave packets moving in opposite direc- 
tions and winding over the circle. We express gQ^{t) as a series in winding number 
and derive tractable approximations for each term. This allows to give qualitative 
explanation and quantitative description to various finite-size effects such as partial 
revivals and regular-to-chaotic transition. 
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1 Introduction 

Exactly solvable spin chains are widely used as toy models for exploring various aspects 
of quantum dynamics. Recent progress in experimental techniques allows to construct 
quantum systems with effective spin chain Hamiltonians (see e.g. [1]), which opens 
new prospects for exploring fundamental concepts such as thermalization and emergence 
of chaos, as well as for applications such as quantum state transfer through quantum 
wires [2j. This motivates further efforts to understand dynamics of spin chains in detail. 

We consider the reduced dynamics of a single spin in the cyclic spin 1/2 XY chain with 
finite number of spins, N. Initially one selected spin has a given polarization while other 
(A — 1) spins are completely uncorrelated and unpolarized. We study the 2;-component 
of polarization of a spin as a function of time. It can be expressed through the two- 
spin time-dependent correlation functions g^{t). Although many papers starting from 
the pioneering paper on XY chain [3] were devoted to calculation of various correlation 
functions, most of the studies concentrated on the thermodynamic limit A — )■ oo. Exact 
expression for g^{t) in the XY model with finite N was derived in IH [3]. It involves 
sums of ~ A oscillating terms and thus is not easily tractable. However, it allows to plot 
g^{t) for various values of model parameters. Such plots readily reveal a rich variety 
of spin evolution patterns which call for explanation (see figures in the present paper, 
especially fig. [T|. One striking feature of the evolution is regular-to-chaotic transition: 
g^{t) is described fairly well by A — )■ oo approximation (which is given by a rather regular 
function of time for a wide range of model parameters) up to some threshold time tth? but 
at tth this concordance is abruptly destroyed by sharp revival; at later times the evolution 
becomes less and less regular and ends up with apparently chaotic fiuctuations near the 
long-time average. This feature is common for all finite spin chains; in particular, it was 
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observed in the XX (isotropic XY) model |6l Ej, in the XX Z model with long-range [8] 
and nearest-neighbor [9] couplings, in the XY model [10]. 

What underlies the regular-to-chaotic transition in the cyclic chain is the winding of 
two oppositely directed wave packets created by the initial perturbation [7] (or reflection 
from the ends of the chain in case of linear chain [6]). Threshold time corresponds to 
the time necessary for the forefronts of the wave packets to make one round trip over the 
circle. The interference between the forefronts of the wave packets and their own tails 
produces partial revivals at t = tth, 2tth; ••• and leads to the regular-to-chaotic transition. 

In order to study spin dynamics in finite spin chains at times greater than tth it 
is desirable to have tractable approximate expressions for g^{t) valid for t > tth- The 
goal of the present paper is to obtain such approximations for the cyclic XY model. The 
mathematical method which we introduce is closely related to the physical picture of wave 
packet winding over the circle and in fact allows to quantitatively describe such winding. 
A winding number s given, we are able to get the approximate expression which is valid 
up to (s -|- l)tth and describes interference between the components of the wave packets 
which have completed 0, 1, 2, s round trips over the circle. These approximations are 
fairly accurate even when the evolution is already apparently chaotic. A related result 
was obtained in [11] in a special case of XX model: they represent Green function of a 
quasi-particle through a sum over winding numbers. 

The obtained approximate expressions involve ~ s oscillating terms and therefore are 
much more tractable than the exact formula as long as s <^ A^. They allow to look at the 
regular-to-chaotic transition (as well as on some other peculiar features of spin evolution 
in finite chains) from a new perspective. However in the present paper we only briefiy 
discuss this issue postponing more extended study for future work. 

We also touch the issue of incomplete thermalization of spins in the XY spin chain. In 
particular we show that the autocorrelation function go^{t) at infinite temperature never 
changes its sign in contrast to what should be expected in case of complete thermalization. 
This intriguing property was previously proven in the special case of the XX chain [6l [7] 
and observed in numerical calculations of spin evolution in the XXZ model with long- 
range interactions [H]. 

The rest of the paper is organized as follows. In Sec. [2] we briefiy describe the XY 
model on a circle. In Sec. |3]we discuss exact formula for gn^{t) and rewrite it in a form 
suitable for subsequent analysis. In two special cases (one of them is the XX model) 
we represent g^^ (t) in transparent and convenient way through the infinite sum of Bessel 
functions; such representation allows to obtain simple successive approximations valid up 
to times tth, 2tth, •••• In Sec. |4] we present our main result - the successive asymptotic 
approximations in a general case. In Sec. |5] we discuss the emergence of chaos in the 
course of spin evolution. The results are summarized in Sec. [6] Bulk of technical details 
is presented in Appendices. In Appendix |X] we describe the diagonalization of the XY 
model. In Appendix |B] we rederive the exact formula for g^{t) at infinite temperature 
using a method which is somewhat more direct than one implemented in the original 
work |5] . These two appendices mostly contain widely known calculations and results; 
however we include them in order to introduce our notations, to emphasize some salient 
features usually omitted in the literature and for the sake of completeness. In Appendix [C] 
the dependance of the wave packet forefront velocity on the model parameters is inves- 
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tigated. Appendix [D] contains the technical details of calculating asymptotic expressions 
presented in Sec. |4} 



2 XY model on a circle 

We consider a chain of coupled spins 1/2 with the following Hamiltonian [21 



N , N 



H = I Y.({i + 7)«+i + (1 - ix^i^i) + ^ E <• (1) 

n=l 71=1 

Here the index + 1 is identified with 1, and is supposed to be even. Two parameters 
enter the Hamiltonian, the anisotropy parameter 7 and the magnetic field h. Without 
loss of generality one may assume 7 > 0, /i > (see Appendix |A]) . In Sec. |4] we will 
concentrate on the case > 1, 7 G [0, 1]. 

Finite XY model is "almost diagonalizable" through the sequential Jordan- Wigner, 
Fourier and Bogolyubov transformations [31 IH IS]- The details of the XY spin chain 
diagonalization are presented in the appendix [Aj 

An important property of the XY Hamiltonian is that it commutes with the parity 
operator H = Y[n=i "^^^ diagonalization leads to the "almost free fermion form" of 
the Hamiltonian: 

H = p^^^ Yl - ^) + p""^ E - (2) 

where 

_ Ar N N ^_r^l^3 N 1 

Aodd - {- Y + 1, + 2, yl, Xev -{-— + -,-— + -,...,— - -\, (3) 

{cg, q G Xodd} and {Cg, q G Xev} are two sets of fermion operators (note, however, that 
two operators from different sets do not satisfy fermion anticommutation relations, see 



eq.(A-21)), p°'^'^ and P^^ are parity projectors, 

P<=" = (1 + H) /2, P°'^'i = (1 - n) /2, (4) 
and fermion energy is defined as Eg = E{(p{q)), ip{q) = 

E{ip) = ^/e{^y^^^Til^, e{ip) = h - cosip, r(v9) = 7 sin v?. (5) 

One can see that the Hilbert space is divided into a two subspaces with odd and even 
numbers of fermions correspondingly. Number of fermions is an integral of motion, and 
when it is fixed, the model looks like a free-fermion model. 



3 Reduced dynamics of a spin at T = oo 
3.1 Correlation function: sum over modes 

We focus our study on the z-component of the n'th spin polarization vector as a function 
of time: 

p:(t)=tr[p(t)<], (6) 
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Figure 1: Patterns of spin dynamics for various values of model parameters. Exact auto- 
correlation function gQ^{t) is plotted. According to eq. ^ it is equal to the polarization 
of the first spin pl{t) provided pl{0) = 1. Horizontal line (blue online) marks the long-time 



average of the autocorrelation function given by eq.(B-lO). Number of spins here and in 
all other plots in the paper is = 100. 



where p{t) = e *'^*p(0)e*^* is the density matrix of the whole chain. 
We choose the following initial condition: 

p(0) = 2-^(li + pi(O)o-i) ® l23...^. (7) 

It describes a situation when at t = the first spin has an arbitrary polarization pi(0), 
while other [N — 1) spins are completely unpolarized and uncorrelated. If the first spin 
is regarded as an open system, while other {N — 1) spins - as an environment, then such 
initial condition corresponds to infinite temperature of the environment. Given the above 
initial condition, the polarization p^{t) can be expressed through the two-spin correlation 
functions at infinite temperature, p^(t) = Pi{0)g^'^i{t), where 

gT{t)^2-''tTK^,{t)al]. (8) 

Due to conservation of parity g^{t) = g^{t) = and we are left with 

p:it)=pm9n-iit). (9) 

Note that such relation between the polarization of a single spin and the correlation 
function holds only in the case of infinite temperature. 
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Thus our problem reduces to investigation of the zz correlation function. Due to 
integrability of the model it may be calculated exactly [H E]. For the completeness of the 
presentation we provide the details of calculation in Appendix [B] The result reads: 



jn 2 /^n 

' '-"odd 



2 2\ 



where 



COS nip{q) cos Eqt, 
^ cos rnp{q) sin Egt, 
sin rnp{q) sin Egt. 



(10) 



^ev{odd) (^) 


= A^- 


' E 






-^rv{odd) (^) 


= A^- 






cjeX(,v(odd) 


^ev{odd) (^) 


= A^- 


E 




cjeX(,v(odd) 



(11) 



In what follows we mainly concentrate on the evolution of the first spin which is 
distinguished by the initial condition. It is described by the autocorrelation function 

As was noticed in [7j, in case of the XX model (7 = 0) Qnit) is always non-negative 
(because C'CT(odd)(^) ~ ^) other words, spin polarization never changes its sign. We 

see that this is not the case for an arbitrary site n in a general XY chain. However, 
the polarization of the first spin still never changes its sign since C'ev(odd)(^) ~ ^ f*-*^ 
7. Intriguingly, the same property (non-negativity of go^{t) at infinite temperature) was 
observed in numerical simulations for the XX Z model with long-range interactions |H]. 
This suggests that this effect could be generic for a large class of spin systems. 

Surprisingly enough, evolution of spin polarization described by the exact formula ( |lO| 
exhibits a rich variety of patterns depending on h and 7. Examples are given in Fig. [Tj We 
aim at explaining major features of evolution and at providing tractable approximation 



to eq.(lO). 



3.2 Correlation function: sum over winding numbers 



Let us rewrite formulae (11) for n = in a different form: 



where 



^°,,(t) = A,it) + 2ZMt), Alit) = Ao(t) + 2E(-l)^AW, 

i=i i=i 

00 00 

= Boit) + 2j:B,{t), Blit) = Bo{t) + 2j:{-iyB,{t), 



A, it) 



TT 

(27r)~^Re j e 



i{E{v)t-jNip) 



(12) 



(13) 



— IT 

As will be shown below j corresponds to a number of windings of a forefront of a wave 
packet produced by the initial perturbation. To obtain the above expressions one should 
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Figure 2: Qo^it) for the XX model. Threshold time is tth = 100. Solid line corresponds 
to the exact expression. Dotted line (magenta online) corresponds to the approximation 
(17) with j = 0,1, ...,10. One can see that the approximation starts to deviate from 
the exact expression only at t ~ lltth- Approximation obtained in thermodynamic limit 
(j = 0) is also shown (red online). It accurately describes gQ^{t) up to the threshold time. 
Horizontal line (blue online) marks the long-time average of the autocorrelation function. 



take a discrete Fourier transform of the r.h.s. of eq.(ll) and use 



1 if / = jN, J G Z, 
otherwise, 
(-1)^- if / = jN, J G 
otherwise. 



(14) 



Formulae (12) have an important advantage compared to formulae (11): infinite sums 



in (|12j) may be truncated at some small j to obtain excellent approximations for times 
t < {j + l)tth with threshold time tth ~ N. Thus one may deal with only few terms in eq. 



12h in contrast to N terms in (11). This statement will be proved in full generality in 



what follows (see Sec. |4]and especially Appendix D.2.2). In two special cases described 
below one can check it immediately. 



3.3 Special case: XX chain 

When 7 = functions Aj and Bj can be expressed through Bessel functions of the first 
kind: 

Aj{t) = (-l)^J'/2cos(/lt)J,jv(t), 



(15) 



JjN{t) is negligible for t < jN, which justifies the truncation of the sums in (12). Threshold 
time in this case equals A^. 



In fact in the case of XX chain eq. ( 10 ) may be further simplified to obtain 




jjN{t)jj'N{t)- 

j+j'=0(mod2) 



(16) 
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Note that h drops out from the final expression. This can be easily seen from the defini- 
tion (pi of Qn^it) if one recalls that | ^„cr^ commutes with the total Hamiltonian 



Eq. (16) may be used to obtain successive approximations: 



'0 



it) 



+ 4J2 



N 



it) 



+ 4J2 



2N 



{t)+4Mt)J2N{t) 



for t e 



[0,N) 
[N, 2N) 
[2N, 3N) 



(17) 



The first line (j = 0) represents a well-known result obtained in thermodynamic (A^ — oo) 
limit [12]. Approximations in which (s-|-l) Bessel functions are kept correspond to s round 
trips of a spin wave over the circle. We postpone further discussion of the physical sense 
of the obtained results to the next section. Exact and approximate expressions for gl^{t) 
in the XX chain are plotted at Fig. [2] 

Closely related results for the XX model were obtained in [TT] . Namely, they represent 
one-particle Green function, which is in fact equal to the zero temperature correlation 
function Qn^if) |t=o = (ii ••• i l'^n+i(^)'^i'l a ■■■ \), as an infinite sum of Bessel functions. 



3.4 Special case: Ising chain with h = 1 

In the case /i = 1, 7 = 1 one obtains 

^j{t) = <^2iAf(2t), ^^g^ 
^jit) = 2 ('^2jAr+l(2t) — J2jAr-l(2t)) = — J2j^(2t), 

where prime stands for the derivative. Again Aj{t) and Bj(t) are negligible for t < jtth 
with tth = N. Successive approximations can be written analogously to the XX case 
discussed above. 



4 Asymptotic approximations 

In the present section we derive asymptotic approximations for functions Aj,Bj which 



enter eq.(12). As we will see, these approximations physically correspond to taking into 
account spin waves which wind over the circle j times. The details of the calculations are 
presented in Appendix |D} Here we outline only major results emphasizing their physical 
meaning. In the present section we restrict our study to the case > 1, 7 G [0, 1]. 



4.1 Winding of a wave packet over a circle 

We approximately calculate Aj(t) and Bj(t) using the method of the steepest descent in 
the plane of complex variable (f. The saddle points for Aj are obtained from the equation 

v{^)t-jN = 0, (19) 

^As Prof. Perk noted in private communication, another way to explain this fact is to use the trans- 
formation into the rotating frame, a procedure famiUar in the theory of magnetic resonance. 
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values of 


7 = 


7 = 


= 1 


h = 


h = 


1 


h=l, 


parameters 




h < 1 


h > 1 




e [0, 3/4) 


1' e [3/4, 1] 


7 = ^2 - 


cos ipo 





h 


1/h 




272+1-1/472+1 


1 


~ 0.414 


~\/l+7 


2(1-72) 


V 


1 


h 


1 


1-7 




7 


2(V2-1) 



* - bulky (although explicit) expression. 

Table 1: Value of maximal group velocity V in some special cases 



where v{ip) = d^pE is the group velocity corresponding to momentum (the equation 



corresponding to Bj is slightly different, see eq.(D-3) in the Appendix D). Note that in 
general the positions of saddle points depend on time (to be more exact, on the ratio 
t/jN). Two important cases should be distinguished, t < jtth and t > jtth, where 
tth = N/V and V = supv{ip) = v{ipo). In the former case eq.(19) has no real roots and 



as a consequence Aj{t) and Bj{t) are severely suppressed (in accordance with a general 
result [13] )• This explains why one can keep only j terms in eq.(12) whenever t < jtth- If 
h is not too close to 1, the suppression law reads 



Aj{t), Bjit) ~ exp [-const ■ jN 



Jt 



th 



t 



Jt 



th 



t < jtth. 



(20) 



where the constant is of order of one and depends on h and 7, see Appendix |D.2.2 

In the opposite case t > jtth eq.(19) has two real roots and Aj{t), Bj{t) are not 
suppressed. 

Threshold time tth is a time which is necessary for the fastest spin wave to make 
one round trip over the circle [7]. Thus Aj(t) and Bj(t) describe contributions of those 
parts of the wave packet which have completed exactly j round trips over the circle. 
The propagation of the wave packet is visualized in Fig. [s] (see also an analogous figure 
for the XX model in [7]). As was shown in j7], initial perturbation localized in the 
first spin gives rise to two wave packets which travel in opposite directions. Each wave 
packet is a superposition of all spin waves of corresponding direction. The velocity of the 
forefronts of these wave packets coincides with the maximal group velocity of the spin 
waves V. Therefore as long as t < tth = N/V, the wave packets propagate as if the chain 
were infinite, and the evolution of the first spin is described merely by oscillations in the 
common tail of the wave packets. This stage of evolution is the only one which may be 
catched by the N ^ 00 approximation. Mathematically it is described by keeping only 
j = terms in eq.(12). 

At t = tth the forefronts of two wave packets complete the round trip over the circle 
and meet at the first site. At this moment the regular evolution of the polarization of the 
first spin is abruptly interrupted by a partial revival. The succeeding evolution between 
tth and 2tth is determined by the interference between the fastest parts of wave packets 
which have already made one round trip and the common tail of the wave packets with 
zero velocity which still stays at the first site. Mathematically this stage is described by 
keeping j 



0, 1 terms in eq.(12). 
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Figure 3: Propagation of two oppositely directed wave packets along the spin chain. 
Snapshots of polarizations of all spins at times t = 5,25, ...185 are presented. In this 
figure the initial perturbation is localized at the 51'st site (in contrast to the rest of the 
paper), the zeroth site is identified with the hundredth one. The model parameters are 
= 1, 7 = a/'v^I'-^, the threshold time is tth — 117.7. Note the emergence of flat 
polarization plateau which grows up to tth/2, then shrinks and completely disappears 
at threshold time. This feature is specific for the case h = 1. A three-dimensional plot 
representing a more generic case (without plateau) of propagation of wave packets along 
the cyclic XX chain may be found in [7j. 
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Subsequent stages are described in a similar fashion. The wave packets continue to 
wind over the circle. At stth < t < {s + l)tth polarization at the first site is a result of inter- 
ference of waves which completed 0, 1, s round trips over the circle. This corresponds 



to keeping j = 0, 1, .., s terms in eq.(12). The revivals at t = stth become less pronounced 
with increasing s due to the decrease of the maximal amplitude and the smearing of the 
forefront of the wave packet, see Fig. [3] 

Clearly the maximal group velocity V is an important quantity in the above picture as 
it determines the threshold time. We show in Appendix [C| that V G [2{^/2 — 1), 1] as long 
as one restricts himself to the case h > 1. Without this restriction V is confined to the 



interval [0, 1]. V = 2(a/2-1) is achieved at/i = l, 7=Va/2-1 (this case is presented on 
the upper right plot in Fig. [T]). More detailed considerations including several important 
specific cases may be found in Appendix [C| 

The above described physical picture of propagation of wave packets and emergence 
of revivals implies that the forefront of the wave packet is rather sharp. This is indeed 
true for the following simple reason. As long as cpo is a point of maximum, a bunch of 
fermions exists with (f{q) lying in the vicinity of <^o- The group velocities of these modes 
are equal to each other and to the maximal velocity V up to quadratic terms. It it is 
exactly these modes which form the sharp forefront of the wave packet which smears very 
slowly compared to the rest of the wave packet. Some proposals for high-quality quantum 
state transfer along spin chains exploit this feature (see e.g. [14j). 



4.2 Asymptotic approximations for t < tth- 



First we separately consider the case t < tth (see Appendix D.l). This case is special 
because the position of saddle points do not depend on time. Approximations at this 
time interval coincide with formulas obtained in thermodynamic limit. 

Let us introduce a dimensionless parameter e = h — 1. Here and in what follows we 
mainly concentrate on the case where h is not too close to 1. In this case the method 
of the steepest descend can be applied straightforwardly, and we obtain an asymptotic 
approximation for times max{l, e~^} < t < tth '■ 



1 TT 

go"{t) - ^ ■ ((ao+ - ao-f + 4ao+ao- cos\t - -)), (21) 



where 



h±l 



/i± (1 -72 



Note that e should be greater than N ^, otherwise the time interval at which the approx- 
imation is valid vanishes. This restriction is relaxed if 7 <^ e. In the latter case even for 



e ^ 1 formula (21) is valid for I <^ t < tth- In fact in this case one may approximate 



(y'g^(t) simply by the autocorrelation function for 7 = which is given by Jo(^) f*^^ ^ < ''-th 



according to (17). 

In the case of small e the application of the method of the steepest descend is more 
sophisticated. The complications are not unexpected because h = 1 is the point of 
quantum phase transition (QPT) for the XY model. However given certain relations 
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Figure 4: Autocorrelation function g^{t) at t < tth or, equivalently, in thermodynamic 
limit. Solid line - exact expression, dashed line (magenta online) - asymptotic approx- 
imation (21) valid for lager times, e"^ < t < tth; dotted line (red onhne) - asymptotic 



approximation (22) valid for smaller times, 1 <t <^ 76 



between e, 7 and ^ one may still obtain accurate approximations, see Appendix D.l 



For example in the case < 7^ ^ 1 we are able to obtain an asymptotic expression valid 
for 1 < t < 7^e"2 : 



1 



2-7^ 



7r(2 -72)t 



2v^ 



■ exp[ 



27^ 



+2, 



2-72 



exp[ 



2x7(1 -7') v^r^ 



-t\ cos(2t — -7 — arctan 



r 



(22) 



If 7 is not too small, the exponents rapidly decrease with time and one is left with non- 
oscillating decay: 



1 



7r(2-72)t' 



(23) 



For certain values of parameters both asymptotic approximations (21) and (22) may 
be applicable, but at different time intervals. An example is given in Fig. |4j 
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Figure 5: Exact autocorrelation function gQ^{t) (solid line) and the approximation cor- 
responding to 5 complete round trips over the circle (points). One can see that the 
approximation excellently describes both the revival and the apparently chaotic evolution 
far from the revival. 



4.3 Asymptotic approximations for t > jtth- 

Now let us turn to asymptotics for functions Aj{t) and Bj{t) and corresponding ap- 
proximations for gQ^{t) in case of j > 1. It was already noted that Aj{t) and Bj{t) are 
suppressed for t < jtth- 

For sufficiently large times and h not too close to 1 (more explicitly, for {t — jtt^) ^ 1 



and e ^ <^ jtth) we obtain (see Appendix D.2.1) 

1 




cos{tE{ipi) — jNifii + 



4' 



(24) 



+ 



\E"{^2)\ 




B,{t) 




\E"{^2)\E{ip2, 



sm.{tE{ipi 



?,m{tE{if2) - iNip2 - 



jN^i + -) 



(25) 



Saddle points v^i,2(^) are obtained from Eq. (19). The latter may be reduced to a poly 



nomial equation (D-16) of fourth degree with regard to cosy). 



These asymptotics being plugged into eq.(12) excellently approximate gQ^{t) every- 
where but in the vicinity of points jtth where revivals occur. In order to describe revival 
one should use different approximate expressions presented in the next subsection. 

The case of small e is again more cumbersome. However it can also be treated as is 
demonstrated in Appendix |D. 2.1 
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4.4 Partial revivals 



Above derived approximations based on the method of the steepest descent are not apph- 
cable in the vicinity of multiples of tth when two saddle points are close to each other and 
to the point cpo. However, as long as at t = jtth is exactly the time when a partial revival 
occurs, it is highly desirable to have an approximation which works well for t ~ jtth- We 
present such an approximation in the Appendix |D.2.3[ It is based on the fact that in 



the case under consideration the integrals in the definitions (13) of Aj and B- pick up 



the major contribution in the vicinity of (fQ. This justifies the expansion of E{(p) in the 
vicinity of ipo which leads to the desired approximate expressions. If h is not too close to 
1, namely e ^ 'fN~^, they read 



\E"'{'Po)\t 



E{<fio) \\E"'{<fio)\t 



■ Ai 



■ Ai 



-N 



tth \\E"'{,po)\t 



ttb \\E"'{^o)\t 



(26) 



where Ai(x) is the Airy function of the first kind and ifo corresponds to the maximal 
group velocity. Curiously enough, the above approximation works well even far from jtth- 
Note that as long as E"'{ipQ) does not depend on time in contrast to E"{(pi^2), it is much 
easier in practice to calculate the r.h.s. of equations (26) than the r.h.s. of eqs. (24), 



(25). The only disadvantage of the approximation (26) is that we do not analytically 



control the errors of this approximation; however, numerical calculations show that they 
are small. 

To summarize, in order to approximate the autocorrelation function up to (s + l)tth, 
one should take Aq, Bq according to eq. (D-7), Aj, Bj with j = 1, 2, s — 1 according to 



eqs. (24), (25) and Ag, Bg according to eq. (26). The resulting expression approximates 



the autocorrelation function with excellent precision as shown in Fig. |5} 

The case when /i ~ 1 is as usual more cumbersome. We do not provide a complete 
analysis which would be rather bulky, however we derive an approximation for = 1, 7^ > 



3/4, see Appendix D.2.3, eqs. (D-32)-(D-36) 



Let us discuss the law which governs the decrease of revival amplitudes. In general 
the s'th partial revival is described by the the mutual interference between all Aj{t) with 
j = 0, l,...,s and mutual interference between all Bj{t) with j = 0, l,...,s. However 
as a first approximation one may consider only As{t) and Bsjt) which give the leading 
contribution. Aj{t) in eq.(26) decreases as t~^^^, while in eq. d24| - as This means 

that the amplitude of revivals decreases more slowly than the averaged value of go^{t) 
between revivals. In fact this makes the revivals so visible against the background. As a 
result for h — 1 > 'j/N one gets from eq.(26) the following law: 



9o 



'th revival 



[sN) 



-2/3 



(27) 



This law work satisfactory for sufficiently large number of spins and for moderate s. 
In particular, as long as the long-time average of Qq^ is of order of N~^, this law can 
not be valid for s > V^- In fact it breaks down somewhat earlier because at large s 
contributions from Aj, Bj with j < s start to contribute significantly. Our numerical 
calculation show that for 10000 spins the above law is reliable for a few dozens of revivals. 
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Figure 6: Extremely pronounced revivals occur in otherwise apparently chaotic regime 
in the small region of the parameter space in the vicinity of the point /i = 1,7^ = 3/4. 
Threshold time here is tth — 117. 



For more moderate number of spins, ~ 100, the law is quickly distorted due to the 
above mentioned contribution from Aj{t) and Bj{t) with j < s. In particular, maxima of 
revivals do not decrease monotonically in this case, see Fig. [6j 

Noteworthy, when h = 1 and 7^ = 3/4 the amplitudes of the revivals decrease even 



more slowly than implied by eq. (27), namely 



n 

i/O U'th revival 



{sN)-'/',h = l, 7^ = 3/4, 



(28) 



see Appendix D.2.3 and especially eq. (D-36) for the details. Again this law works for 



sufficiently large number of spins. However the fact that the h = 1, 7=3/4 point of a 
parameter space is a special one reveals itself already for modest N ~ 100: the revivals 
appear to be especially pronounced in the vicinity of this point (see fig. IgI), although 



they do not decrease monotonically due to the above discussed interference of Aj, B 



with different j. A very similar effect was observed previously in [15]. Namely, it was 
numerically discovered that when a perturbation propagates from one end of an open- 
ended XY chain to another, the attenuation of the amplitude of a wave packet is minimal 
for /i = 1, 7 ~ 0.7. 



5 Transition from regular to chaotic evolution 

Plots of g^{t) presented in the present paper clearly demonstrate that transition from 
regular to chaotic evolution is a general feature of spin dynamics. In the present section 
we provide a discussion of this fact on a qualitative level. More thorough study including 
quantitative considerations will be presented elsewhere. 

From Fig. |2]it is evident that the spin evolution is apparently regular at small times but 
apparently chaotic at large times, threshold time tth determining the relevant timescale. 
However it is not so easy to define the terms "regular" and "chaotic" rigorously in the 
present context. One should especially be cautious when using the term "chaos" here. It 
is well known that the notion of chaos in quantum systems is not uniquely defined (see 
e.g. [16]). One widely used definition of quantum chaos is based on energy level repulsion. 
According to this particular definition XY model is certainly not chaotic because it is 
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Figure 7: Sensitivity of the correlation function to the small variation of 7. Function 
plotted is the absolute value of difference of two correlation functions gQ^it) correspond- 
ing to two slightly different values of anisotropy parameter: 1(70^(^)^=0.5 ~ 5'o^(^)7=o.505| • 
Magnetic field in both cases is h = 1, number of spins = 100. Threshold time is approx- 
imately 112. The horizontal line (blue online) shows the value which the above difference 
would admit if the two functions were absolutely uncorrelated. 



integrable and thus its level statistics is Poissonian (i.e. non-repulsive). In what follows 
we briefly discuss two distinct approaches which may be used to describe the level of 
irregularity of Qnif)- 

The first approach is based on the physical picture of winding of the wave packet over 
the circle and exploits asymptotic approximations derived above. In this approach we 
pragmatically consider the evolution to be regular at some interval of time if the correlation 
function may be well approximated by a linear combination of few (^ A^) oscillating 
functions with different frequencies (probably multiplied by a power-law prefactor) at this 
interval. Conversely, the evolution is considered to be chaotic when the approximation 
involves many (~ A^) harmonics. The first stage of evolution {t < tth) is the most regular 
one: according to eqs. (21), (22) it is described by a single cosine. At times tth, 2tth5 3tth 
new functions Aj,Bj come into play in eq.(12) and the number of harmonics increases 
stepwise. Thus the level of irregularity also increases. This does not last forever: according 
to eq.(lO) A^ harmonics is enough to describe go^it) exactly. Evidently the largest possible 
level of irregularity is achieved not later than at t = A^tth- Note that here we use the term 
"harmonics" in a slightly non-standard way: we do not demand that the corresponding 
frequencies should be multiples of a single, minimal frequency. The described approach 
resembles the Feigenbaum rout to chaos through period doubling (see e.g. [IZj). However 
in the case under consideration there is no doubling - the relation between frequencies of 
new harmonics switching on at certain times is not so evident. Moreover, these frequencies 
may even be slowly varying in time. The Fourier analysis of the correlation function at 
different time intervals is necessary to obtain more quantitative picture. This will be done 
elsewhere. 
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In the second approach one examines the level of sensitivity of g^{t) to small variations 
of Hamiltonian parameters h and 7. This approach introduced in [18] resembles the 
definition of classical chaos through the extreme sensitivity to initial conditions. We 
visualize the sensitivity of g^{t) to small variations of 7 in Fig. |7| One can see that during 
the regular stage of evolution (t < tth) such sensitivity is small, while at large times {t > 
few-tth) it is comparable to what one could expect if two correlation functions with slightly 
different parameters were absolutely mutually uncorrelated. As in the previous approach, 
the rate of thus defined chaoticity increases stepwise at times which are multiples of t^h, 
the first step occurring at t = tth being especially pronounced, see Fig. [7j Curiously, 
our numerical experiments indicate that the sensitivity of g^{t) to variations of h and 7 
generically tends to be larger in the vicinity of quantum phase transition line h = 1. It 
would be quite surprising if this relation between QPT and chaoticity of the correlation 
function is confirmed, since the correlation function is calculated at infinite temperature 
while QPT occurs at zero temperature. 



6 Discussion and conclusions 

If one merely plots the evolution of spin polarization in the finite cyclic XY chain for var- 
ious values of Hamiltonian parameters, one immediately discovers that the spin dynamics 
exhibits a variety of peculiar features. Previous studies [HI [7] and the present study reveal 
the following physical picture. 

• A threshold time exists up to which the polarization of a given spin evolves as if the 
chain were infinite. This is the time necessary for the fastest spin wave to make a 
round trip over the cyclic chain. Up to the threshold time the evolution is regular 
(for h>l). 

• At the threshold time the regular evolution is interrupted by a partial revival. Subse- 
quent partial revivals occur at 2tth; Stth? ••• Generically the evolution becomes more 
and more irregular (chaotic) after each partial revival. This picture may be ex- 



plained on the basis of eqs.(12) in which terms Aj{t), Bj{t) start contributing to 
the sums at t = jUh- Physically this picture may be understood as the interference 
between the forefront of the wave packet which have made j round trips over the 
chain and its lagging tail. 

Our main result is as follows. We have expressed the autocorrelation function gQ^{t) 
through infinite sums of functions Aj{t), Bj{t) over the winding number j. This expression 
allows to write down successive approximations for gQ^{t). For any given j the jth approx- 
imation takes into account j complete round trips of the wave packet and consequently 
may be applied up to t = (j + l)tth- Functions Aj(t) and Bj(t) are defined in integral 
form. In two special cases (7 = and 7 = 1, h = 1) they can be expressed through 
Bessel functions. In a general case we are able to provide very accurate approximations 
to these functions valid at various times and in various regions of parameter space. These 
approximation are obtained using methods of complex analysis, in particular, the method 
of the steepest descent. 
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We have shown that the amphtude of the revivals decreases slower than the average 
amplitude of the correlation function between two revivals. This fact makes the revivals 
clearly visible against the background. We have analytically shown that the most pro- 
nounced revivals are observed in a vicinity of a special point of the parameter space 
[h = 1, 7^ = 3/4). Previously this was observed numerically in an open-ended chain [15] . 

A striking feature of the dynamics in the finite spin chain is the transition from regular 
to chaotic evolution. In the present paper we have restricted ourselves by brief and 
qualitative discussion of the nature and origin of this transition. Further work is necessary 
to give a more exhaustive and quantitative analysis. 

We also have observed a feature concerning long-time properties of evolution. It was 
shown in [7j that in the XX model the polarization of any spin never changes its sign 
provided the initial conditions are specified by eq. ([T]). We point out that in the XY model 
this property still holds for the first spin (distinguished by the initial conditions) but does 
not hold for any other spin. Thus the first spin retains the memory of the fact that it was 
distinguished by the initial conditions forever. Moreover, its time-averaged polarization 
differs from the time-averaged polarization of any other spin which was already noted 
in the context of the XX model [7]. Thus we encounter the absence of the complete 
thermalization which is, however, only a finite-size effect (scaling as 1/A^). 
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A Diagonalization of finite cyclic XY spin chain 

A.l Ranges of parameters 

Let us rewrite the Hamiltonian we are going to diagonalize: 

N , N 

H{h, 7, = ^ + + (1 - ^X^n^i) +2^<^ ^^-1) 

n=l n=l 

where indices 1 and + 1 are identified, and A^ is even. Here we have introduced coupling 
constant k, which is taken to be 1 everywhere in the article but this subsection. Let us 
show that one may consider /i, 7, k > without loss of generality. This means that one 
can change the sign of each constant by means of local unitary transformation U. These 
transformations correspond merely to rotations of the coordinate systems at each spin 
site. 
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To change sign of h one can transform a"^ —(^n-> ~^ each spin site n: 

N N 



n=l n=l 

[/t<t/ = <, U^alU = -al, C/t<[/ = U^E{h,^,K)U = 7, /^)- (A-2) 

Analogously, to change sign of 7 one transforms cr^ — >■ cr^, cr^ — >■ — cr^ at each site n by 

AT 

means oilJ — W e^K'^l'^, 

n=l 

To change sign of k one transforms crf^ — > — crf^, crf^ — )■ — crf^ at each even site 2m by 

Af/2 

means of t/ = H e^'^sm'^/^. 

m=l 

As soon as sign of k is unimportant, one may put k — 1. 

A. 2 if in terms of 

We define the operators in a usual way, 

'^^^IK + O, <^n=^«-0- (A-3) 

These operators are neither Bose nor Fermi operators: 

<^i(^n = (^n(^m for m^n. (A-5) 
The following simple equalities prove to be useful: 

= 2(7+(7- - 1 = -2(7-(7+ + 1 (A-6) 

a^a'^ = -a^a^ = cr+, a''a~ = -a" ^ -a~ (A- 7) 

The Hamiltonian may be rewritten in terms of as follows: 

H^Ho + H^ + Hh (A-8) 

with 

1 ^ 



n=l 
N 



Hj = l J2(<j>:!:+i + <yn<+i). (A-io) 

n=l 



n=l 
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A. 3 Jordan-Wigner transformation 

Define the operators 

n 

n=l 

Evidently, IlAr coincides with the parity operator 11 defined in Sec. |2] 
Define Fermi operators a~, as follows 

= O-nHn-l = I^n-lC^n ^ = = ^n-l<^n (^"13) 

This implies 

a~ = a'Un-i = Iln-ia~ , cr+ = a+n„_i = n„_ia+, (A-14) 

{am^an} = Smn, {a+ , 0+} = {o;;, } = 0, (A-15) 

< = 2a+a- - 1 = -2a; a+ + 1. (A-16) 

The Hamiltonian takes the form (note that now the ordering of a^, a^+i is important; 
also note the change of the total sign): 



^0 = -^E(«n«n+i + fln+ifln) " (1 + ^){a%a^ + a|a^)]. (A-17) 

n=l 
N 

^7 = "2 [2^(^n"n+i + «n+i«n) " (1 + n)(a>^ + a^a^)]. (A-18) 

n=l 

TV 

Hh = hY^ a+a- - Nh/2. (A-19) 

n=l 

A. 4 Fourier transformation 

Define for arbitrary real q 

h- = ^—y e-2-^("-^)/^a;, 6+ = Y e2-'?("-i)/^a+, (A-20) 

* n=l * n=l 

1 1 _ p27ri{fc-g) 

{&^^n = K-,&,} = 0, {^^,^-} = ^ ^_^,..(._,)/^ - (A-21) 



Then 



In particular, if one takes 



or 



g = -Y + l, -y + 2,..,- (Xodd) (A-22) 

N I A^ 3 AT 1 , 

^ = -y + 2' -y + 2'-'y-2 ^^-^ ^^-'^^ 
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then the set of is the set of Fermi annihilation operators. 
The Hamiltonian may be written in terms of 6^ as follows: 

H = i/°ddpodd ^ ^evpev ^^,34^ 

with 

podd ^ (1 _ n)/2, = (1 + n) /2, (A-25) 
N/2-1/2 N/2-1 

i^-= Yl ^<?' E H, + Ho,N/2, (A-26) 

9=1/2 g=l 

H, = {h- cos<^(g)) (6+6- + + 7sin<^(g) (6+6^^ + 6:^6") - /i, (^(?) = 27rg/7V 

(A-27) 

and 

^o,iv/2 = (/i - l)hth^ + {h + - /i. (A-28) 

A. 5 Bogolyubov transformation 

Define the following quantities 

r, = 7 sin ip{q) , = /i - cos ip{q) , = + (A-29) 
Each i7q may be written as follows: 

It is possible to diagonalize this matrix through Bogolyubov transformation 
The diagonalization condition reads tan^^ = ^q/^qi we choose 



Cq = cos ^ h- + sin ^ ht,. (A-31) 



r 

= arctan ^ for all q^Q. (A-32) 

This transformation preserves the anticommutation relations. i?o,Ar/2 requires special 
treatment, which leads to 6^/2 = 0, 



0, h>l, 

TT, < /i < 1. 



(A-33) 



= cos ^ - sin ^ clg (A-34) 



The inverse transformation reads 

^q - . ^q 

The odd and even parts of the Hamiltonian take the form 

^odd(ev)^ E E.ic'q^--^). (A-35) 

(ev) 

This completes the diagonalization. 
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A. 6 Eigenstates 

Let us first prove the existence of tlie Fock vacuum states witli respect to tlie anniliilation 
operators c~, i.e. tlie states |vac)odd5 |vac)ev wliich satisfy 

Cg |vac)odd (ev) = V g e Xodd (ev)- (A-36) 
Evidently it is sufficient to prove that 

n ^ 0- (A-37) 

"jG^odd (cv) 

If this condition is fulfilled, one can always choose some states |\l'odd(ev)) and normalization 
constants Kev such that 

|vac)odd = ^oddC_^y2+l'^-Af/2+2'"^Af/2l^odd); (A-38) 

|vac)ev = ^cvC_jy^2+l/2'^-Af/2+3/2'"'^A'/2-l/2l^ev)- (A-39) 



The equality 



{'^-7V/2+l' ['^-Af/2+2' {'^Af/2-l' ['^7V/2 Y\ '^q^^'''^^~^ (A 40) 



and the analogous equality for q E Xodd (ev) prove eq.(A-37). Note that |vac)ev is indeed 
an eigenstate of the Hamiltonian, while |vac)odd is not. 

All the eigenstates of the Hamiltonian are obtained from the vacuum states by applying 
the creation operators c^. To create the odd number of fermions one should use q G Xodd 
and |vac)odd, while to create the even number of fermions one should use q G and 
|vac)ov 

Evidently one can enumerate all the eigenstates of the Hamiltonian by the multiindexes 

Qm = {qi, q2, Qm}, 0<M<N, (A-41) 
with the ordering qi < q2 < ... < qM- Then an eigenstate with M fermions reads 

\Qm) = c+ ...C+C+ |vac)odd(ev) (A-42) 



with qi,q2, ...,qM £ Xodd (ev) when M is odd (even). The corresponding eigenenergy reads 

^Qm E (A-43) 

(?eQi\/ '?e-'^odd (cv) 

For our purposes we need only the matrix elements between the states with the same 
parity, therefore we use the notation |vac) without subscripts in what follows. 
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B Calculation of g'^if) 

To calculate the correlation function at infinite temperature, 

9::it) = 2-^5^(Q|a^|Q)(g|<^JQ)e-^(^'^-^'5)*, (B-1) 

one needs to calculate the corresponding matrix elements. To do this one uses 

COS I COS ^ctcp e2^^(P-P)"/^ + sin | sin ^c^c+e-^^^^P-P)"/^. (B-2) 



Here p, p can run either through Xodd or through - the expression is valid in both 
cases. Now it can be easily seen that only three types of matrix elements do not vanish: 

1. Diagonal matrix elements. 

(QK+ilQ) = ^ E ^(Qm.p) cose,. (B-3) 

Here ri{QM,p) = 1 if p G Qm and —1 otherwise; Xm = -^odd (ev) if M is odd (even). 

2. Matrix elements between two states with equal number of fermions, differing by one 
fermion momentum. 



Qm = Km-i U {p}, Qm = Km-i U {p}, p,p ^ Km-i, p p : 

2 9r, -\- 9. 
N 



{QK+M = e^-(^-W^^ COS ^^(g|ctc,-|Q), (B-4) 



where {Q\cpCp \Q) = ±1, depending on the signature of the corresponding permu- 
tation. Note that this sign is not important for calculation of g^^it). 

3. Matrix elements between two states one of which can be obtained from another by 
addition of two fermions. 

Qm = Qm-2^ {p} ^ {p}, p,p^Qm-2, P^P- 
(QK+ilQ) = (QK+ilQ)* = -|e^-^(^+^W^sin^^(Q|crc;|Q). (B-5) 



Now let us sum in eq. (B-1) separately over each type of matrix element: 



1. Sum over Q = Q gives 



2-N 



\ P,peXodd odd M Qm / 
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2. Sum over pairs (Q, Q) of the form Q = K U {p}, Q = K U {p} gives 



E + E 

C V f 



,2ni{p-p)n/N-iiEp-Ep)t ^^^2 ^P + % 



3. Sum over pairs {Q, Q) of the form Q = K U {p} U {p}, Q = K gives 
4Ar2 



27ri(p+p)n/Ar-j(£p+£;p)t ^^^2 ^P ^ 



(B-7) 



(B-8) 



while summation over Q = K, Q = K U {p} U {p} gives a complex conjugated 
contribution. 



If one takes p = g in expression (B-7), it becomes equal to the (B-6) contribution. 
Exploiting this one readily obtains 



9n{t)= E + E 

vP.gGXodd p,q&x 



cos 



COS 



f27r{p-p)n \ 2^P + ^P 

( {Ep-Ep)t ) cos 

2tt{p + p)n 



N 



{Ep + Ep)t sin 



It can be straightforwardly verified that this expression leads to eqs. (10), (11). 



Eq. (B-9) can be used to find a long-time average of the autocorrelation function 

T 

Qq^ = lim J gQ^{t)dt. Let us assume that there are no degeneracies in Eg other than 







the mirror degeneracy Eg = E^q (in other words, that Eq = Ep implies |g| = \p\). This is 
a generic case. Then 



N 2N 



(B-10) 



VpSXodd psx 

e 



Term —1/N in parenthesis emerges due to g = 0,N/2. Long-time average of correlation 
function for 7^ can be calculated analogously. In the specific case 7 = the result 



(B-10) coincides with the expression obtained in 



C Group velocity of spin waves 

In the present section we consider /i > 0, 7 G [0, 1]. Group velocity of spin waves reads 
v{ip;h,'j) = {h — (1 — 7^) cos v?) sin Lp/E{Lp; h, 7). (C-1) 
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We are interested mainly in the maximal velocity for given values of parameters h and 7 



V{h, 7) = supv{<^; h, 7) = v{<^o; h, 7), (C-2) 



where ipo is the supremum point. Due to the symmetry of E{ip) we can consider > 
without loss of generality. Extremum condition |<^o ~ leads to the fourth degree 
polynomial equation 

P{z) = {l-byz^-3h{l-b)z^ + {2b{l-b) + h'^{3-2b))z'^-h{h'^ + b)z-b{l-b) + h% = 

(C-3) 

with z = cos !fQ and 6 = 7^. We are interested in the real roots of this equation which lie in 
the interval [—1, 1]. Let us show that there is only one such root whenever h > 1 (this fact 
is important for the application of the method of the steepest descent, see Appendix [D|. 
In this case the above equation implies that z > 0, therefore in fact we have to consider 
the interval [0,1]. Since -P(O) > 0, -P(l) < and P(+oo) = +00, we could have 1,2 or 
3 roots in [0, 1]. If there were 2 or 3 roots of P = in the considered interval, then the 
equation P' = would have 2 roots in [0,1]. However the latter equation has no more 



than 1 root in the considered interval (z= ^^-^_^^ ). Thus equation (|C-3|) has exactly 1 root 
in the interval [0; 1] for h > 1. 

Let us now consider several important special cases. 

I. 7 = 0. In this case ipo = 7i/2, V = 1. 

II. 7 = 1. In this case 

( h, h < 1 A\ 

cosy^o = \ u-i A ^ 1 (C-4) 

and 

y = <! ? : t ; (c-s) 



7 





h < 1 


I h- 


\ h>l, 


h, 


h < 1 


1, 


h>l, 


V = 


1-7. 



III. h = 0. In this case cosy^o = ■ y 

IV. h = 1. This is an especially interesting case as it corresponds to the quantum phase 



transition. Velocity v{(p) has a step at = 0, step height being equal to 27. Eq. (C-3) is 
simplified to 

{z - 1)2 ((6 - 1)2^2 + (252 _ 5 _ 1)^ + 52^) ^ Q_ 
One should distinguish two cases. 

IV a. 6 G [3/4, 1]. In this case the only root that satisfy \z\ < 1 is z = 1. Thus ipo = 0, 
r = 7. 

IV b. 6 G [0, 3/4). In this case cosy^o = '^^'^2(i-b)^ from which one can easily write down 
an expression for V which appears to be somewhat bulky. One can also find a minimal 
value of V with respect to 7 : 



inf \/(l, 7) = 1/(1, V2 - V2) = 2(v^- 1) = 0.828427... (C-7) 

7e[o,i] " 

In what follows we show that this is the minimal value of V in the whole region /i > 1, 7 G 
[0,1]. 

Y. h ^ 00. In this case y^o — ?■ 0, F — )■ 1. 
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Let us investigate how V{h, 7) varies with h. The derivative over h has a rather simple 
form: 

dhV{h, 7) = b{l-h cos ¥?o) sin yDo/E^{<^o; h, 7). (C-8) 

To calculate it we used that dhV{h, 7) = dhv{(po; h, 7) due to the equation d^v{ipo; h, 7) = 0. 
The stationary points of V{h,'y) with respect to h are given by dhV 
cos(y9o = 1/h. We plug the latter equality into eq. (C-3) and obtain 



which leads to 



(1-7^)(1-.W-(1-7V)=0. 



(C-9) 



The only roots that satisfy < 1 are z = ±1 which correspond to h = 1. This point 
is not extremal because V{0,j) < ¥{1,^) < \^(+oo,7). Thus for any fixed 7 maximal 
group velocity V{h, 7) monotonically grows with h from 1 — 7 at = to 1 as — )• 00. 
As a consequence, if one considers only h > 1, than the minimal value of V is given by 
eq. 



D Asymptotic expressions 

Here we consider in detail asymptotic expressions for spectral functions Aj{t) and Bj{t). 
Let us explore domains in which E{ip) is an univalent analytical function. The branch 
points of this function are found from the equation [h — cos(p)^ + (7sin<^)^ = 0. For 
/i^ > 1 — 7^ its solutions are ipbr± and (plr±, where 




cpbr± = 27r/c + i ■ arccosh , k E Z. (D-1) 

\ 1-72 

The domain where E{(p) remains univalent analytical function is the hole complex plane 
without row of branch cuts from ip^r- to ipbr+ in the lm{ip) > half-plane and row of 
branch cuts from ipl^_ to ipl^_^^ in the lm{{p) < half-plane. It is important that lm{E{Lp)) 
is positive in I and III quadrants and negative in II and IV quadrants. Functions E[(p), 



e{(p), exp[—ijN(p], as well as integrands in eq.(13), are 27r periodic functions of Re{(p). 
Thus we consider the strip Re{{p) G (— tt; tt]. 



In order to use the method of the steepest descent for the integrals in eq.(13) we have 
to find saddle points for functions 

fAj{v>) = iE{v)t - tjN<f, fBj{^) ^ iE{^)t - ijN^ + In (D-2) 

The saddle points are defined by the equations 

E't-jiV = 0, and iE't - ijN + {ln^ = (D-3) 

correspondingly. 
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D.l Asymptotics for spectral functions of zero order 

Let us first consider spectral functions of zero order, Agit) and Boit). There are four 
saddle points in the strip Re{ip) G (— 7r;7r] which for Aq read 

ifi = 0, = TT, Lfs = —i ■ arccosh -, (^4 = +i ■ arccosh -. (D-4) 

1—7^ 1 — 7^ 

For Bq we have exactly the same points yjf = y^i, = ip2 and slightly (for t ^ 1) shifted 
points and (pf: 

f3A = =F« ■ arccosh (^^3^ + ^3,4 ■ ^ + (^-5) 



One can find 03^4 by substituting ip^^ = 993,4+ A in the second equation in (D-3). Note 
that if in the (c3,4-t~^)-vicinity of 993,4 the series for E{ip) is convergent, then the difference 
between 993 and ipf^ is of order of t~^. 
Let us define the following parameter: 

e = h-l. (D-6) 

For large enough e we can find asymptotics in the case t 3> max{e~^, 1} in a straightfor- 
ward way. Indeed, we can transform the integration path from our initial Cq (integration 
along Re((/9) = from (p = — vr to ip = ir) to the path Ci which goes through I and III 
quadrant, where lm{E) > 0, and through saddle points (pi = and ip2 = n (see figure 
left). Then we immediately have 



Ao{t) - ^2^(ao_cos((/i-l)t + f) + ao+cos((/i + l)t-f) + 0(f)), 
Bo{t) - ./^(ao-sin((/i-l)t+f) + ao+sin((/i + l)t-f) + 0(i)). 



(D-7) 



where 



h±l 



/i± (1 -72) 



and we take into account that e{0)E^^{0) = 6{'k)E^^{7i) = 1. Region of applicability for 
large enough e is given by t ^ 1, and for e — )■ we have t > 2e~^ In ^ where A stays for 
the value of an error. Under these conditions we have the following approximation for 
t < tth: 

1 TT 

9o"it) - ■ ((^0+ - ao-f + 4ao+ao- cos2(t - -)). (D-8) 
Znt 4 

Numerical evolution shows an excellent coincidence with the exact solution in the region 

1 ^ t < tth- Note that this expression becomes asymptotic for Jo(t) in the case 7 = 



(XX chain) which is in accordance with eq. (17). 



Let us consider the case of such small e that e ^ 7 and e <^ (t)~^- Since we are 
interested in the dynamics on time scale of order of tth or larger, the latter condition in 
fact implies e <^ A^~^ ^ 1. Now we cannot integrate over the contour Ci due to the 
small convergence radius of the series in the vicinity of 991 = (note however that the 
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Figure 8: Integration paths for Aq, Bq in the case t 1, t ^ e^^ (left) and in the case 
^ t ^ 1, 7 ^ e (right). Empty triangles correspond to saddle points, filled circles 
- to branch points of E{ip). 

contribution from the saddle point = n remains intact). Therefore we use another 



integration path, C2 (see fig. (8j)), which goes through saddle points — vr, ips, ipi, ip^ and 
ends up in ip)2 = tt. Consider integration path from ip^ to ^94. Integration over branch cuts 
does not contribute to Aj, Bj (this statement is true for spectral functions of all orders). 
This is because Iie{E{ip)) = on the branch cuts, and therefore 

Re [ e'^^^^^-'^'^dip = lm [ ^M^e''^^^'>-'"^dip = 0. (D-9) 

icut J cut E{ip) 

For small e branch points can be expanded as ipbr-.vlr- — ±^-e7~^+0(e^7~^). At the seg- 
ment v^br-] functions i?((/9) ande(v9) are real- valued, and £'(</)) < -\/2e, e{'p>)E~^{Lp) = 
1 + 0(e), therefore 

fPbr- _ 

Re e'*^(^)rf¥P < (727-^6 + 0(eV'))-min{et,l}, (D-10) 

Ini ^e'm^)d^ = 7-ie + 0(e^-2). 

Combining these results with contributions from saddle points ip2, fs, ^i, one obtains 



;i _ exp|--i=i] + 0(i) + 0(67-') I . 



For Bq, keeping in mind (D-5), one obtains 
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+(1 - f )-i exph— 2L=t] + O(^) + 0(67"^^ 



Region of applicability for these asymptotics is limited by the condition that exp[—^\E"{ip3) 
must be small enough (here R is the radius of convergence for E{(f) near ips^^). We have 
R = |v?3 — ^br+\ — (2 — \/2)7 for small 7. Thus we get the following condition of applica- 
bihty for the above asymptotic: 

t>(2-y2)-V'ln(^), (D-13) 

where A is the order of the relativ e err or. U nder this condition we can neglect the 
term exp[— 7^(1 — 7^)~2t] in (D-11) and (D-12). Numerical evaluation shows excellent 



coincidence of these asymptotics with the exact values of Aq and Bq in the case 7 -C 
1, t ^ 1. Moreover, these expressions exactly coincide with asymptotic forms for Bessel 



functions in (15) in the case 7 = 0, h = 1 which is not obvious from our derivation 



method. With these remarques, we find 



2-72 272 
g'^U) = — — 1 + — -^=exp[ 1] 



/ 2 - 7^ -y2 ^ ^ 

+2. / — , exp[ t1 cos(2t arctan ) 1 , (D-14) 

which gives us an excellent approximation for QQ^if) in the case e <^ 7, e ^ A^~^, t < tth- 
For not very small 7^ we can neglect exponential suppressed terms and obtain 

D.2 Asymptotics for spectral functions of non-zero order 

For spectral functions of non-zero orders the situation is slightly more complicated. Again 
there are four saddle points in the strip Re((/9) G (— tt, tt], but now their positions vary 
with time (see fig. [9]). Consider the function Aj^ j > 1. The corresponding saddle points 
satisfy the 4-degree polynomial equation on 2; = cos ip : 



(1 - J^Z^ - 2h{l - 72)^3 + [((/i2 _ (1 _ ^2)2)) ^ _ y)]^2^ 

+ [h{l - 72) - 2Ch]z -e + C{h^ + 7^) = 0, 



where ( = {Nj)'^{t) . This equation viewed as the equation on ip gives eight solutions 
in the strip Re(9?) G (— 7r,7r]. Four of them are relevant (i.e. are solutions of eq. ( |D-3 )) 



and other four are irrelevant (are solutions of the equation —E't — jN = 0). Since E'((p) 
takes all possible real values on the each branch cut, one pair of saddle points, ip^ and 
ip4, lies on two brunch cuts symmetrically with respect to the real axis, analogous to the 
j = case. 

Let us consider the positions of two other saddle points, ipi and ip2, especially their 



evolution with time. The definition of the threshold time implies that for t < jtth eq.(D 



has no real roots. Thus ipi and ip2 are complex. In fact they are complex conjugate 
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Figure 9: Integration paths for Aj, Bj in the case t > jtth (left) and in the case t < jtth 
(right). Empty triangles correspond to saddle points, filled circles - to branch points and 
to the point ipo which satisfies E{ipo)" = 0. 



to each other. When time goes on they both approach (po which lies on the real axis, 
and eventually merge ai t = jtth '■ V'i(j^th) = ^2{jtth) = V'o- For t > jtth fi and 
lie on the real axis and move apart from ipo and from each other, approaching and vr 
correspondingly as t — )■ oo. 



D.2.1 Asymptotics for t > jtth 

To start with, we warn the reader that we do not provide a strict mathematical proof 
that the suitable integration path exists which goes through the chosen saddle points 
in all presented cases. However, the existing of these paths looks quite natural in all 
cases and, moreover, corresponding asymptotic expressions show excellent coincidence 
with numerical evolutions. Strict mathematical proof is postponed for further work. 

With this warning made, let us turn to the case t > jtth- We start from the case 
of non-small e and we assume that ipi and !f2 are situated far enough from each other, 
so that we can neglect their mutual influence in the asymptotics. The conditions under 
which this assumption is fulfilled are considered in what follows. 

Note that we can integrate along the path Ci U Co U C2, where Co is the original path, 
Ci starts from = —n — ioo and goes to = —it along Re*/? = const = — tt, C2 starts from 
Lf = TT and goes to (p = n — ioo along Rey? = const = vr (see fig. [o]). Since fA,Bj{x + iy) are 
2tt periodic functions of x, the value of the integral along the new path is exactly the same 
as along Cq. Now we transform the path Ci U Co U C2 to the path C which starts from 
Lp = — TT — ioo, goes through saddle points (pi and ip2 and ends at (p = n — ioo (see fig. [9} 
left). There are two topological different possible integration paths C: the first one goes 
above the branch cut, while the second one - under the cut. In the latter case according 
to the Cauchy theorem we have to subtract the integral over the branch cut. However as 



was discussed above this integral does not contribute to Aj and Bj (see eq.(D-9)). 



Now we can proceed to find asymptotic expressions as contribution from points ipi 



and ip2. Under all specified conditions, we immediately obtain eq. (24) for Aj. For Bj 
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using reasonings similar to those following eq. (D-5), we get eq. (25). 



Let us now investigate the range of applicability of eqs. (24), (25). Firstly we consider 



in what cases we can use the standard approximation for contribution of saddle points 
under the assumption that radius of convergence for corresponding series is large enough. 
In this case the derived approximation may deviate from the exact expression for two 
reasons: small value of \E{(pi^2)"\ and interception of contributions for (pi and ip2 due to 
their close relative position. These two features can appear only for small times after jtth- 
Let us give more precise estimation without detail explanations. If 6t = t — jtti^ > 0, then 



it has to be 

tth 



> 



where r = 



is quantity of order of 1 for the vast majority 



of the Hamiltonian parameter space. One can see that these asymptotic approximations 
for spectral functions of order j become accurate starting from time close to jtth- 

Now let us investigate in what cases series does not converge in large enough circle for 
some saddle point. We have to explore small enough e, at least ^ t. Since for spectral 
function of order j we are interested in t > jt^h, it is useful to consider ^ jN. Firstly 
we define the position of ipo : 



(po = arccos 



^0 



. 27^ + 1- y/yTT ^ 3 
. 2(W^ ) + ^^^)' ^ ^4' 



(D-17) 



„ 1 1 

2262 



Let us consider the case 7^ < |. For times which satisfy jtth < t < jtth = jNj'^ + 0(e) 



there is no point near Lp = 0, and the derived asymptotic expressions (24), (25) are 
valid. If t > tth one obtains 



^1 



72*2 



1\ PN^ 



- 1 



i>jtth=J- + e . 
7 



(D-18) 



Thus for t ^ e^^ we cannot use the above derived approximations because v?i is situated 
close to yjbr- and the radius of convergence R = \/ e^^"^ + tpl is very small, \E" [Lpi)\tR^ < 
1, thus we cannot use the method of the steepest descent for (pi (we assume here 7 e). 
Instead we can proceed analogously to the case of spectral functions of zero order. Namely, 
we move the integration path in order to go through saddle points ip^ and </94 and neglect 
the value of integral between (plr_ and (pbr- as we have done in (D-7). The difference from 



the case of Aq and Bq is that now we can neglect exponentially suppressed contribution 
from (pz, but contribution from ip^ may be not small for some portion of time. Numerical 
evaluation shows that for ^ 1 the contribution from (p^ is exponentially suppressed at 
a timescale ~ tth- Summarizing, for time t > we obtain 



A,(t) 




E"icp,)\ 



exp[-ti-iE{ip^ + 0)) + jN{-iipi)] 



(D-19) 



+ 



\E"{V2)\ 



COs(tE(v92) - jN(p2 



4' 
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where one should remind that —iE{(p4 + 0) > 0, 
that we have to take only one half of contribution from (^4. Analogously one obtains 



iip4,) > 0. Factor ^ in due to the fact 



^ \ 2 (-ZE((^4 + 0)) V \E"iip4)\ 



exp[-t(-zE((^4 + 0)) + jN{-i^^)] (D-20) 




^sin(tE(^.)-,iV^.-^) + 0(i: 

The only case we do not investigate here is 7 < e <C 1. We leave it for future work. 
If 7^ > I, then ipo is situated in the vicinity of ip = 0, which leads to tth = t 



th 



A^7~^ + 0(e). Therefore asymptotic expressions (D-19), (D-20) are valid starting from 
time t for which <^2 is situated far enough from (f = (the corresponding condition reads 
^\E" [ip2)\t(p2 ^ 1). For times in the vicinity of jtth another approximation should be 



used, see Sec. D.2.3 



D.2.2 Asymptotics for t < jt 



th 



When t < jtth the path of integration should not go through the I quadrant, since for 
all values of s G (— 7r;7r) dyRe{iEt — imp) > (remind that ip = x + iy). We shift our 
integration path to C or C (see figure ^) and obtain asymptotics from the contribution 
from only one saddle point, ipi. For large enough \6t\, where 6t = t — jtth < 0, we have 



2n\E"{ipi)\t 



cos{E{(pi)t - jNipi + (r{ipi)), 



(D-21) 



where (T{ipi) is the angle between the path and the line Im/^ = const. Let us evaluate 
ipi approximately for t < jtth such that t^^6t <^ 1. If Eq = E"'{ipo) 7^ 0, then 



(pi = (po - if] + 



IE, 



(IV) 



6 K 



5 .E, 



{IV) 





- i-X^^)rf + 0{r^^) + iO{r^% (D-22) 



72' E ' 
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where 



2N5t 



\Eq \ttth 



6t 

-2r — , r = 



_K_ 
IK' 



Up to the first order in 5t ■ t ^ c:^ 6t ■ j ^t., one obtains 



A, it) 



ri(jN)^ -5t, 



247r2 



exp[ 



-2^2? 3 
-^—jN{ )-2 



th 



+ 0(( 



th 



-St. 3, 

M 



'Jt 



th 



)]cos{E{^o)t~jN^o) 

(D-23) 



jtth ' 



This expressions gives only the order of suppression; if one is interesting in more precise 
expression, he has to directly solve eq.(D-3) to find an exact value of ipi and substitute 
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it in the general formula (D-21). This formula (and approximation ( |D-23 )) is valid until 
we can neglect the term ^E"'{ipi){A (p)^ in comparison with ^E"{ipi){A series 
expansion near Lpx. This leads to the condition jN{^)^ ^ 1. In the opposite case, 

3 J I J I 

) are not valid and the leading order contribution is 
rich leads to 



< 1, eqs. (|D-3|) and (|D-23 



given by the term |£""((/)i)(A ^p)'^, w 



r(i)r3 



23367rj3 A^s 



mpoY'it 



th 



\EipT)\t 



cos{E{ipi)t - jNpi 



(D-24) 



where values of E{(pi) and E{ip'l') may be found according to eq. (D-22). Eqs. (D-24) 



and (D-23) become asymptotics for Jj{t) in the case of XX chain. If one wants to derive 
asymptotics valid in the region jiV(-^)2 ~ 1, he has to calculate the integral through 



saddle point with proper path direction (here we use approximation (D-22) for saddle 
points), 



' saddle 



exp[--7]\E'^'\tz' - -\E'^'\tz']dz 



and use the formula 



(D-25) 



A,(t) 



27t' 



saddle 



(t) exp[ 



-^—jN{ — )-^]cos{E{po)t - jNpo] 



(D-26) 



Eqs. (D-23) and (D-24) can be obtained from the above formula by neglecting the second 
summand in the exponent in (D-25) and by expansion of the integrand in eq.(D-25) in 
powers of rj. 

Let us turn to Bj. In order to describe its behavior in analogous way one should start 
from 



Bj{t) ~ Im 



27r\E"{ip,)\tE{ip, 



■ exp[iE{(pi)t — ijNcpi + ia{(pi)] 



(D-27) 



instead of (D-21). We do not describe in details Bj{t) because there is no simple ap- 
proximation formula for (D-27) for all possible values of cos{(po). However the exponential 
suppression for Bj{t) has the same form as for Aj{t), only the preexponential factor 
differs. This is because the suppression is determined by the exponent of the quantity 
—Re{—iE{ipi)t + ijNipi) which is the same for Bj and Aj up to t~^. 

To conclude this subsection, the suppression of spectral functions Aj{t) and Bj{t) with 
6t = t — jtth < with exponential precision reads 



Aj{t) ~ Bj{t) ~ exp[ 



2E'{ip, 



JN{- 



"''-)i + o((^; 
jtth 



(D-28) 



This is in accordance with a general result [Kj. Here we assume that E"'{(Pq) ^ 0; for very 
small E"'{lpq) one has to take into account ^''^^^^(v^o) which leads to a slightly different 
law. 
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D.2.3 Asymptotics for t ~ jtth 

The considerations of the present subsection are less rigorous than in the previous ones. 
The only consequence, however, is that we do not completely control errors for derived ap- 
proximations. Numerical calculations demonstrate that the latter are nevertheless rather 
accurate in a wide range of model parameters. For certain regions of the parameter space 
in which the derived approximations fail we are able to identify the reason and point out 
the way to overcome the difficulties. 

Let us describe the method. When t ~ tth the saddle points are situated near ipQ. 
Integrand is a very fast oscillating function for all except ~ (^jq- Thus the value of the 
integral is picked up on a small segment [ipo— A ip; ipo+ A (f], A ip > \{y:>2 — V^i)- In order 
to estimate errors for this approximation, one has to make bulky calculations in the spirit 
of the above subsections. We avoid this in the present work. 

In order to calculate the integral along the small segment, we expand E{(f) in the 
vicinity of <^o- And the last approximation is to replace the interval of integration from 
[— A if, A if] to [—00; 00], where the integration variable is /i = (/? — (/?o- The latter trick is 
justified because our new integrand oscillate as ~ exp(i ■ const ■ fi^) when fi — )■ ±00. All 
these approximations are legitimate when the model parameters are such that ipo is far 
enough from points of branching. Thus the approximation is valid for large enough e and 
all values of 7 and for e <^ 1 it is valid for 7^ < |. Let us assume that Eq = E"'{ipQ) is 
not very small. In this case one can consider power expansion of E{(p) up to {if — ifo)^ 
and neglect terms ~ 0{{ip — ipoY). Thus for t ~ jtth one gets 



^jit) = TT / co^{Et-jN^)d^ ^ — / co^{Eot-jN^Q + {E'^t-jN)^i--\E'^'\t^?)dn 



1 



3 1/" f — if^u / 2 \ 3 1 

/ cos[iEot-jNipo)-—^N(——] ^+-eH- (D-29) 



l-^o'lV 27r tth vl-Eo'lv 3 

Here all functions / with subindex should be understood as /o = f{fo)- The above 
integral is a well-known Airy function of the first kind: 

Ai(x) = - / cos{—+xOd^ (D-30) 



Let us investigate now the case = 1, 7^ > ^ (we do not investigate here the case of 



Thus we get eq. (26) (with analogously reasonings for Bj). 

3 

4 

small, but non-zero e = h — 1 <^ 1). For h = 1, 7^ > | we cannot use formula (D-29), 

since ipo = ipbr- = ^1^- = 0. But we can argue that the main contribution for Aj is picked 

up on the segment [0, A y?], since only on this segment oscillation frequency is not very 

large (note that for [— A ip, 0] the frequency is large; this is due to the discontinuity of the 

group velocity at = 0). If we somehow expand E{ip) on [0, A we can obtain a good 

approximation. Consider -E{(p) on a complex plain. If we have values for E{ip) on the 

segment ip G [0,7r] fixed, we can make analytical continuation to R.e{ip) G [— vr, 0] in two 

different ways: with the branch cut: Re(v?branch cut) = 0, Im(v9branch cut) e [-ipl^^,ipbr+] 

or with the branch cut: Re(v?branch cut) = 0, Im(v5branch cut) e [-00, -ipl^^] U [ipbr+, 00]. 

In the first case we have our original function E{ip) in the segment ip G [— 7r,0]. In the 
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second case we have some new function E{ip) which does not coincide with E{(f) in this 
segment. But in the latter case we can use power expansion for E{ip) in the circle with 
radius R = \(Pbr+\'- 



3 - 3 ^ 167^ - 15 



247 



I92O73 



(D-31) 



Analogously to ( |D-29[ ) one obtains 



1 

2^ 



3-472 

cos[(7t - jN)ip H -tif" + 

247 ^ I92O73 



3 167^ - 15 , 



(D-32) 



If one can neglect the term ~ ip^ (i.e. if 7^ — | is large enough), he gets 



27 



2 \{l'-l)t 



Ai 



-ilt-jN) 



27 



, /i = 1, 7^ > -. (D-33) 



For Bj one uses power expansion e{ip)E ^(i/?) in the vicinity of (/? = to obtain 



27 



47 \{l^-\)t 



■ Ai' 



-ilt-jN) 



27 



(7^ - l)t 



/i = 1, 7^ > -, (D-34) 



where prime stands for the derivative. 



Let us consider the special case h = 1, 7^ = 1. We have tth = ■ The term ~ ip 



gAi„(a;) 



2 



in power expansion E{!^) is zero: -^(v^) = — j^^j^V^^- Let us introduce the function 



TT 



^71 

cos( VxC}di 

n 



(D-35) 



Airy function of the first kind is a particular case of this function: Ai(x) = gAi3(x). 
gAi„(x) for n > 3 exhibits the same behavior as Ai(x): for positive x it is exponentially 
decreasing, and for negative x it oscillates and goes to zero when x —00. Aj and Bj 
are expressed through this function and its derivative: 



A, it) 



1 / 2AV3\ 
-2[—) 



gAis 



VSi /24y3V 



h = l,^' = - (D-36) 



3 
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Note that rather small value of the coefficient of term Lp^ in power expansion of E{i~p) 
implies that these approximations work well for sufficiently large time {t ^ 24\/3). Anal- 



ogously we require t ^ 27(7^ — |) ^ in eqs. (D-33), (D-34). Note that if one wishes to 



investigate approximations which successfully describe Aj, Bj near threshold time in the 
region of parameter space /i = 1,7^ ~ |, then one has to calculate integral (D-32) saving 
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both terms ~ (p^ and ~ (p^ . Therefore, in this case there is no such clear power law for 



time dependence of maximum value of spectral function as in (D-37). 

Now let us discuss possible values of maximums for Aj and Bj. If jN is sufficiently 
large the positions of global maximums of Aj and Bj coincide with those for Ai(x) and 
gAi5(a;) : 



sup Aj it) 



■ 03, for h> 1] or = 1, 7^ < 



(D-37) 



supAj(t) ~ - 
t ^ 



27 



(7^ - t; 



as, for /i = 1, 7^ > 



as, for /i = 1, 7^ 



3 
4 



Here as and as are global maximums of Ai(x) and gAi5(a;) (as = 0, 54..., as = 0.44...) cor- 
respondingly. Without loss of precision one can replace t by j'tth in the above expressions. 
We do not present the analogous expressions for Bjii) because positions of maximums of 
these functions obviously do not coincide with those for Aj (when Ajif) achieves global 
maximum Bj(fy becomes zero). In the last two cases in eq. (D-37) it is Al- which gives 
the main contribution in g^^ii) near t = jtth- When h > 1 or h = 1,7^ < | one has to 
investigate maximum of Aj{ty + Bjlt)"^ in order to find the leading contribution to gQ^{t). 
This leads to 



MaX(revival j) [Qq'' 



MaX(i.evival 
MaX(revivalj)[5'( 



it)] 
it)] 



27" 
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l|)a3 for h > I; or /i = 1,7^ < 
for /i = 1, 7^ > 
for /i = 1, 7^ = 



3 

4' 

3 
4' 

3 
4" 



(D-38) 

Note that this is a rather rude approximation, since interference terms may be large. 
Therefore these expressions work well only for large number of spins. Numerical evolution 
shows that the revivals are maximally pronounced for /i = 1 and 7^ slightly less then |. 



This may be expected on the basis of eq. (D-32). 



Let us emphasize once more time that all the derived expressions can be derived 
with more rigor using integrals in the complex plain analogously to that was done in the 
previous subsections. We have not explored the whole parameter space. In particular, we 
have not considered the cases /;. = l,7^~|,e = /i — l^l,e7^0, or/i> 1,Eq' = 0. 
Some of the derived approximations work well only for large time and, correspondingly, 
large number of spins (for example, t ^ 24-\/3 or t ^ 27(7^ — |)~^). However, for a 
large region of parameter space these approximations work fairly well, and they provide 
an opportunity to investigate amplitudes of maximums in partial revivals, or at least 
the law of there decrease. For these reasons we decided to include in the paper these 
not completely rigorous calculations. The validity of formulae derived in the present 



subsection is justified by the fact that approximations (26) give the same result as the 



more rigorously derived eqs. (D-23) and (D-24) when the ranges of applicability overlap. 
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